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ABSTRACT 

Using the theoretical formalism of C. Regal etal. PRL (2004), we 
have studied the BCS-BEC crossover regime. We have evaluated 
the value of energy gap parameter and chemical potential for the 
crossover as a function of dimensionless parameters IK F aj" 1 from 
BCS-theory. Our theoretical result indicates the crossover occurs 
in the very small region of the parameter -1< IKpaJ" 1 <1. 

Keywords: BCS-BEC crossover, Feshbach resonance, Pseudo- 
gap, chemical potential, Energy gap parameter, evaporative 
cooling. 



1. INTRODUCTION 

In earlier paper 1 , we have studied 
the BCS-BEC Crossover physics and 
evaluated energy gap parameter A/E F and 
chemical potential |o/E F from BCS theory. In 
this paper, we have evaluated the Super fluid 
Density (p s /p), Quasi Particle Contribution 
(p F n /p) and Fluctuation Contribution (p B n /p) 
as a Function of (T/T c ) for BCS limit [(K F 
a s ) 1 = -2.0], pseudo gap limit [(K F a s )" = 
0.0] and BEC limit [(K F a s ) _1 = 2.0]. Here K F 
is Fermi wave vector and a s is the s-wave 
scattering length. 



As we know that in a Fermi gas with 
a Feshbach resonance (FR), one can tune the 
strength of the paring interaction by 
adjusting the threshold energy of FR. 2 The 
BCS-BEC crossover has been realized by 
using this unique property. 3 ' 4 Here, if one 
increases the strength of the paring 
interaction the character of super fluidity 
continuously changes from weak coupling 
BCS type to strong coupling BEC type of 
tightly bound cooper pairs. 5 ' 6 In the super 
fluid phase, the super fluid density p s is the 
most fundamental quantities. The value of p s 
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is always equal to the total carrier density p 
at T = 0 while it vanishes at the super fluid 
phase transition T c . These properties are 
satisfied in both Fermi and Bose super 
fluids, irrespective of the strength of the 
pairing interactions 7 . But there is a crucial 
difference between p s in a Fermi super fluid 
and that in a Bose super fluid. In a mean 
field BCS theory p s (T>0) originates from 
the thermal dissociation of cooper pairs. The 
resulting normal fluid density p„ = p - p s is 
determined by quasi particle excitations. On 
the other hand p n in the Bose super fluid is 
dominated by Bogoliubov collective 
excitation. 7 Therefore it very interesting 
problem to see as to how p s in a Fermi super 
fluid changes into p s in a Bose super fluid in 
BCS-BEC crossover. 



Here one assumes two atomic hyperfine 
states described by pseudo-spin o =t,|. 
¥ p + = (C pT + , Cpj) is a Nambau field operator. 
C po + is the creation operator of a Fermi atom 
and Tj are the Pauli matrices ( j =1,2,3) which act 
on the particle - hole space. 
^ p =e p - (j. =( p 2 /2m -p. ) is the atomic 
kinetic energy measured from the chemical 
potential u, U is the tunable pairing 
interaction associated with FR. 

P P 2 P ~2 (2) 

Pj q represents the generalized density 
operator. p iq and p 2q describe the amplitude 



In this paper taking theoretical 
formalism of Y. Ohashi 8 and Y. Ohashi & 
A. Griffin 9 ' 10 , we have theoretically 
evaluated the super fluid density in the BCS- 
BEC crossover. These authors have taken an 
uniform super fluid Fermi gas at finite 
temperature and extended the strong 
coupling Gaussian fluctuation theory for T c 
developed by Nozieres and Schmitt-Rink 
(N.S.R.) 4 ' 5 to super fluid phase below T c . 
They self-consistently determined A and \i. 
We have used their formalism to calculate 
(p s /p), (p F n /p) and (p B n /p) for BCS-BEC 
crossover. 

2. MATHEMATICAL FORMULAE 
USED IN THE EVALUATION 

One starts with the BCS 
Hamiltonian in Nambau representation 



(1) 

and phase fluctuations of the order 
parameter A respectively. In equation (1) 
the interaction is described by the sum of the 
interaction between amplitude fluctuations 
(Piq Pi -q) an d the phase fluctuations (p 2q p2- q >. 

In the NSR theory 6 , T c is described 
by the Thouless criterion in the t- matrix 
approximation. The resulting equation for T c 
has the same form as the mean -field BCS 
gap equation with A— K). However, in 
contrast to the weak- coupling BCS theory 
(where [i is equal to the Fermi energy E F ), \i 
remarkably deviates from E F in the BCS - 
BEC crossover regime due to strong 
pairing fluctuations. The NSR theory 



A 2 U 
u p p ^ q 
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includes the strong -coupling effect by 
solving the equation of state within the 
Gaussian fluctuation approximation 6 ' 7 in 
terms of pairing fluctuation. 

Now one extends the NSR theory to 
the super-fluid phase below T c .To calculate 
A, one uses the BCS gap equation 

tanh(^£ F ) 



m „ 



1 



2E 



2£ 



•] 



(3) 



where E p =[ (s p -u) 2 + A 2 ] 1/2 is the 
single particle excitation spectrum . In 
equation (3) one eliminates the well known 
ultraviolet divergence by employing a two- 
body scattering length a s 

4xa r -U 



m 



Pc 1 

i-uY— 

p 



where, a s is two body scattering length. 

To calculate u, one considers 
thermodynamic potential Q. 

Density is given by 



(4) 



the 



P = 



-dQ. 



(5) 



Fluctuation contribution to Q (=d Q) is 
calculated from relevant Feynman diagrams. 
Now summing up these diagrams, one 
obtains total densities 



Ifj^dfl m 2e p 



p/=S(l-f)tanh(^)- 

p " (6) 



p F is the number of Fermi atoms 

in the mean field approximation. In equation 
(5 a) the second term describes the 
fluctuation contribution. Z(q,iC0 m ) is the 

matrix correlation function. IX, is the 

generalized density correlation function. V n 

is the boson Matsubara frequency. (O m is the 

fermion-Matsubara frequency. Super fluid 
density in the BCS-BEC cross-over is 
determined as 



Ps=p-Pn 

p is the total carrier density given by 
P = H 1 + ^ll Tr lhG(p,iG) m )] 



(7) 



(8) 



G(p, iG) m ) = G 0 (p, i(O m ) + G 0 (p, im m (p, iCO m )G 0 (p, iQ) m ) 

(9a) 

where 



G(v,i(0 m ; x =io) m - Se p T 3 + At, 



(9b) 



(5a) 



G Q (p,iCO m )h the matrix single particle 
thermal Green's function. Lis the self- 
energy which involves corrections to G G p n 
is the well known BCS normal fluid density. 
p n is calculated both for boson and fermion. 
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Pi = J Z P* % 0 Tr[Gso(G(p,iCO m )]Q -> 0- G so l = {iCO m - - £ p T 3 + At,) d D 



Q„ 



2m > 



(10) 



3m p 3£ p 



(12) 



Here G 0 is replaced by G s0 .Super current 
state is described 11 by order parameter. 

A(z) = Aexp(iQz) 

Super fluid velocity V s =Q/2m where f(E) is the Fermi . Dirac distribution 

G s is the matrix single particle function. Boson normal density (fluctuation 

Green' s function in the super current state. correction) is given by 



Pn =-| X Pz % Q WGso(pJa)J^(pJa)JGso(p,ia> m m ^ 0 



(13) 



In the Weak-coupling BCS regime 
(K F a s )" 1 « -1, pairing fluctuation are weak 

and one finds that p n _ p F n or p s = p- p F n . 

In this regime equation (12) shows that p n is 

dominated by the quasi-particle excitations 
with excitation gap A . In the BCS - BEC 
crossover regime, the chemical potential 
deviates from the Fermi energy E F and 
becomes negative in the strong coupling 
BEC regime 46 . One can calculate the 
chemical potential jU in the BEC limit 
where 



(K F a s ) »1. Using equation (3) fl is 
calculated as 



M = 



-1 



2ma, 



(14) 



In BEC regime the chemical 
potential /n works as a large expectation 
gap therefore quasi particle excitation as 
well as p F n are suppressed. This shows that 



Cooper pair do not dissociate in the Fermi 
atoms due to large binding energy. 

From equation (13) one can 

calculate the fluctuation contribution p^ . 

This is the dominant term in the strong 
coupling regime BEC. From (13) one 

obtains p% as 



Pn 3M^ q dE B q 



(15) 



where n B (E) is the Bose distribution 
function. M= 2m is the molecular mass 



2m 2m 



(16) 



Equation (16) is the Bogoliubov Phonon 
spectrum in a dilute molecular Bose gas with 
a repulsive interaction V B - 4jr(2a s )/M and 
the BCS order parameter (p-{&J% 7tmA) 1/2 . 
In the BEC regime the normal fluid density 
is dominated by Bogoliubov collective 
excitations in a molecular Bose super-fluid. 
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Table Tl 

An evaluated result of ( P± ) as a function of for BCS limit [(jf F a s )"' =-2], pseudo 



gape[( J ST F a i )- 1 = 0] and BEC limit [(/sT^)" 1 = 2] 









/' 




P 




BCS =- 2 ] 


nseudo eane \(K n V 1 — 01 




O.U 


1.0 


1.0 


1.0 


0.1 


0.975 


0.982 


0.995 


0.2 


0.956 


0.967 


0.977 


0.3 


0.932 


0.955 


0.964 


0.4 


0.897 


0.902 


0.912 


0.5 


0.824 


0.855 


0.866 


0.6 


0.746 


0.797 


0.805 


0.7 


0.618 


0.639 


0.652 


0.8 


0.546 


0.568 


0.574 


0.9 


0.348 


0.382 


0.403 


0.95 


0.226 


0.267 


0.288 


1.00 


0.059 


0.122 


0.147 


1.05 


0.002 


0.097 


0.106 








Table T2 




F y 

vn evaluated results of quasi particle contribution iPn s as a function of r_\ for BCS lim 

P T/ 
(K F a s y l = -2] . pseudo gape^A^)" 1 = 0] and BEC limit [{K F a s y' = 2] 


T 






p.: 




I ) 
P 


BCS [(K F aX =-2] 


pseudo gape [(K F a s ) 1 =0] 


BEC l(K F ay=2] 


0.2 


0.004 


0.002 


0.0 


0.4 


0.098 


0.008 


0.0 


0.6 


0.185 


0.086 


0.0 


0.8 


0.274 


0.105 


0.0 


1.0 


0.456 


0.126 


0.0 


1.2 


0.684 


0.149 


0.0 


1.4 


0.756 


0.185 


0.0 


1.5 


0.889 


0.225 


0.0 


1.6 


0.954 


0.246 


0.0 


1.7 


1.038 


0.278 


0.0 


1.8 


1.176 


0.304 


0.0 
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Table T3 

An evaluated results of quasi particle contribution , P» \ as a function of r 1 \ for BCS limit 

P T c 
[(K F a s y [ =-2],pseudogape[(i: Faj )- 1 = O] and BEC limit [(K F a s y 1 = 2] 



c 


- <"•') - 

P 


BCS [(K F a s r'=-2] 


pseudogape [(K F a s )~' = 0] 


BEC [(K F a s y'=2] 


0.2 


0.002 


0.004 


0.056 


0.4 


0.006 


0.016 


0.089 


0.6 


0.009 


0.108 


0.126 


0.8 


0.012 


0.148 


0.248 


1.0 


0.122 


0.288 


0.336 


1.2 


0.142 


0.326 


0.409 


1.4 


0.167 


0.449 


0.526 


1.5 


0.198 


0.567 


0.678 


1.6 


0.207 


0.659 


0.776 


1.7 


0.226 


0.787 


0.892 


1.8 


0.245 


0.896 


0.967 


2.0 


0.268 


1.052 


1.122 



3. DISCUSSION OF RESULTS 

In this paper, we have presented the 
method of evaluation of ratio of super fluid 

A) 

density and total density P as a function 



T 



gap limit l(K F a s y' = 0]. 



of for BCS limit 1 ^"^ "~ ZJ , Pseudo 

1 and BEC limit 
[(K F a s ) =2] 0ur theoretical results 

A) 

indicate that P is larger in BEC limit and 

T 

smaller in BCS limit as a function of ° . 



(^) 



T 

(^) 



P decreases with T < for all three cases. 
The results are shown in Table Tl. 



In table T2; we have presented the method 
of evaluation of quasi particle contribution 



P as a function of T ' for all the three 
limits. Our theoretical results indicate that 

p F 

{ Pn_ ) 

P is larger for BCS limit and smaller in 
Pseudo gap limit. In table T3; we have 
shown the evaluated results of fluctuation 

Pn 



T 



contribution P as a function of 
for the above three limits. Our evaluated 
results show that fluctuation contribution 

Pi 

( P ) is smaller for BCS limit K*,*,)" 1 = -2] 

and larger for BEC limit K^f",) = 2 l From 
the above calculation, One observes that if 
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one increases the strength of the pairing 
interaction BCS -type normal fluid density 
dominated by quasi - particle excitation 
changes into BEC type normal fluid density 
dominated by Bogoliubov collective 

excitations. As S plays an important role in 
two fluid hydrodynamics these results would 
be useful in study of dynamical properties in 
the BCS-BEC crossover region at finite 
temperature 12 15 . 
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